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Abstract 

We re-examine the well-studied one dimensional spin-1/2 XY model to reveal its nontrivial 
energy spectrum, in particular the energy gap between the ground state and the first excited state. 
In the case of the isotropic XY model - the XX model - the gap behaves very irregularly as a 
function of the system size at a second order transition point. This is in stark contrast to the 
usual power-law decay of the gap and is reminiscent of the similar behavior at the first order phase 
transition in the infinite-range quantum XY model. The gap also shows nontrivial oscillatory 
behavior for the phase transitions in the anisotropic model in the incommensurate phase. We 
observe a close relation between this anomalous behavior of the gap and the correlation functions. 
These results, those for the isotropic case in particular, are important from the viewpoint of 
quantum annealing where the efficiency of computation is strongly affected by the size dependence 
of the energy gap. 
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I. INTRODUCTION 


The ground state of a quantum many-body system changes drastically as parameters in 
the Hamiltonian are driven across a quantum transition point. This phenomenon is preceded 
in a finite-size system by closing of the energy gap between the ground and first excited states 
as a function of the system size growing toward the thermodynamic limit. The rate of the 
gap closing is an important measure characterizing a quantum phase transition. 

According to the finite-size scaling theory, when applied to a second order phase transition 
with a divergent correlation length, physical quantities generally behave polynomially as a 
function of the system size [l-3]. In contrast, the gap is expected to close exponentially 
fast at the first order quantum phase transition: the two ground states at opposite sides 
of the first order transition point have significantly different properties and consequently 
their overlap in a finite-size system is very - typically exponentially - small. The overlap of 
the two states determines the energy gap since the overlap corresponds to the off-diagonal 
elements of the effective two-level Hamiltonian describing the system around the transition 
point and the gap is directly related to the magnitude of the off-diagonal elements. It 
is therefore expected that the order of quantum phase transition in the thermodynamic 
limit is generally in one-to-one correspondence with the rate of the gap closing toward the 
thermodynamic limit, polynomially for the second order transition and exponentially for the 
first order transition [4j, 5]. 

Studies of the energy gap are also important from the viewpoint of quantum annealing 


12 }, in which parameters of the system are controlled in time and driven across a quantum 


critical point. The rate of the gap closing directly affects the efficiency of quantum annealing. 
An exponential closing of the gap implies an exponentially long computation time whereas 


a polynomial gap leads to a polynomial time 
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19]. 


Although the above-mentioned correspondence between the order of quantum phase tran¬ 
sition and the rate of the gap closing generally holds true in most systems, an interesting 
counterexample has been found by Cabrera and Jullien 20, ]2lJ who showed that the first 
order quantum phase transition in the one-dimensional transverse-field Ising model with 
antiperiodic boundary condition is accompanied by a polynomial closing of the energy gap. 
See also Ref. 22] for essentially the same result. Furthermore, another very anomalous 


example has been given for the infinite-range quantum XY model, where the energy gap 
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behaves in many different ways at first order quantum phase transitions |23J. It has been 
shown there that many types of the gap closing: polynomial, exponential, and even factorial 
closing coexist along a critical line in the phase diagram. This example suggests that we 
should be very careful while relating the type of quantum phase transition to the rate of the 
gap closing. 

In the present paper, we first analyze the one-dimensional quantum isotropic XY model 
with s = 1/2, which exhibits he second order phase transition, and show that this model is 
another example of the anomalous gap behavior. Although the properties of one-dimensional 
quantum spin systems have been studied from a number of different perspectives 244371]. 


our focus is on the anomalous system-size dependence of the gap, at variance from most of 
the previous studies. We conclude that the energy gap behaves quite peculiarly (highly oscil¬ 
latory) as a function of t 


quantum XY model 


re system size, in a very similar way to the case of the infinite-range 
231] , This observation should be taken seriously from the perspective 


of efficiency of quantum annealing. 

Next, we analyze the anisotropic XT' model and show that the energy gap also behaves 
anomalously as a function of the system size in some regions of the parameters. Finally, 
we observe that a direct relation holds between the energy gap and the correlation function 
in the ground state. For instance, the anomalous oscillatory behavior of the gap precisely 
coincides with the oscillations of the correlation functions. 

The organization of the paper is as follows. In Sec. II, we define the model, diagonalize 
the Hamiltonian and obtain the energy gap. In Sec. IH. we analyze the size dependence of 
the energy gap using both analytical and numerical methods. In Sec. IV, we show that a 
relation is established between the energy gap and the correlation function in the systems. 
Finally, our conclusion is given in Sec. V. 


II. THE HAMILTONIAN AND ITS SPECTRUM 

We study the one-dimensional quantum XY model with s = 1/2 in transverse and 
longitudinal fields, F and h , respectively, and periodic boundary condition, 

TV TV TV 

H = - 2 + (! - 7K ? <i] - r 5>i - 5>?. (!) 

i= 1 i= 1 i —1 
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where 7 controls the anisotropy, and erf, erf and erf are the Pauli matrices acting at site i. 
For the moment, we assume that the system size N is even. Without loss of generality, we 
assume that Y > 0 and 7 > 0. 

We consider quantum annealing protocol where the longitudinal hied h is driven in time 
for a given, fixed value of the transverse held T. For 7 > 0, there is a hrst order quantum 
transition at h = 0 between the phases with (er f ) > 0 and (erf) < 0 (for h < 0 and h > 0 , 


respectively), as long as 0 < Y <1, see Fig. 1 [271-34]. This expectation value has a jump 



FIG. 1: The phase diagram of the XY model in Eq. (JT]) , which undergoes a hrst order 
phase transition across the surface marked in pink between the phases with positive and 
negative expectation values of erf. Along the isotropic line of 7 = 0 marked in green, for 
0 < T < 1 , the transition is of the second order. 


at h — 0, which is the reason why we call it the hrst order transition. In contrast, for 7 = 0, 
a second order transition takes place, as the expectation value (erf) continuously changes 
across 0 at h = 0 


Our aim is to clarify the system-size dependence of the energy gap at h — 0. Notice that 
the longitudinal held in Eq. (JT]) , — h )TV erf, induces transitions between states with even 
and odd values of the ^-component of the total spin S* ot = f (Tf erf, while other terms in 
the Hamiltonian conserve the parity. For our protocol, it is therefore relevant to study the 
energy gap between the ground and hrst excited states, where one has an odd and the other 
has even value of Sf ot . 

The XY model in Eq. ([1]), for h = 0, can be diagonalized following the standard procedure 
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27H34]. Since the energy gap depends on the subtle hnite-size effects in the energy spectrum, 
we repeat the derivation in some detail. After the Jordan-Wigner transformation, erf = 
2 ejei — 1 , erf — ierf = 2 Q ]f[ ■ <i (—cr|), where <7 are fermionic annihilation operators, the 
Hamiltonian is expressed as 


N—l 


H = 


[( c *Wi - ^ 4 + 1 ) + 7(44+i - c * c *+i) 

Z=1 

+ (-1)^ [(cU - Cjv4) + 7(44 - cjvCi) 

N 1 

-2rJ4(4 Q --), 

Z=1 

where we will call A" c = Y2i c l c i sig 11 operator. After the Fourier transformation, 


c i = 


1 

y/N 




the Hamiltonian in Eq. (j2j) is transformed as 
H = ^H k , 

k 

H k = - \2a\a k {cosk + T) + 7 (e* fc aj.at fe - e~ lk a k a_ k ) 


r, 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


where the values of k depend on the boundary condition: k — k\ for periodic cjv+j = Cj (N ( 
odd) and k = k 2 for antiperiodic cjy + j = —Cj (N c even), with 

2 


2 N 

kl = 0, ±—'. 7T, • • • , ± 71, 7T, 

, 1 3 N — l 

h = ± N*’ ± N*’-' ± —'' W - 


( 6 ) 

(7) 


For 7 = 0, the transformed Hamiltonian H k is already diagonal. The part of the Hamiltonian 
for a given absolute value of momentum reads 


H k + H_ k = — 2 


4 a -k 


cos k + T iy sin k 
—i'y sin k — cos k — T 
We diagonalize it using the Bogoliubov transformation 




l -k 


2 cos k. 


( 8 ) 


1 

^3 


cos y 

• • Qu 

i sin 



1 

_ 1 


i sin y 

Ot. 

cos -f 


a ~k 


cos 9 k = (cos k + T) /e(k), 
sind k = 7 sin k/e(k), 
e(k) = a/ (cos k + T ) 2 + (7 sin k ) 2 . 


( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 


5 























Finally, the diagonalized Hamiltonian in each parity sector reads 
(i) For N c odd 


H °dd = _ 2 ^ r + _ 2 (r - i)4d w + 2 r 



2 53 3 

+ dl fc3 d-fc3 - !)e( fc 3), 

(13) 


k3 




2 4 

N -2 

(14) 

k-i = -7T, -7T, • • • , 

3 N N 

-7T. 

N 

The ground-state energy is 




77’odd 

-^o — 

-2-2^e(fcs)=- 

-2 — 53 e (^’3) + e (0) + e ( 7r ) 



kz 

ki 


j 

\~E kl e ( k i ) 

(r < i) 

(15) 

= 

5 

I 

[ 2 (r-i)-E tl «( fc i) < r>1 ) 


where we have to be careful to pick the state with correct fermionic parity. 


(ii) For N c even 




H mn = -2 y>iA 

+ dl k4 d_ k4 - l)e(fc 4 ), 

(16) 






1 3 

A- 1 

(17) 

Ka = -7T, -7T, 

4 N N 

, -7T. 

’ N 


The ground-state energy is 


= - 2^2 e(fc 4 ) = -£%(**). (18) 

&4 k,2 

The true ground-state energy is given by one of these two possibilities, i?g dd or i?Q Ven , while 
the energy of the first excited state is the other one. Therefore the energy gap is equal 


A(iV,r,7) = |£° dd -£* ven 


k2 ki 

2(F ~ 1) + y; e(k 2 ) - ^2 e (^i) 

&2 k\ 


(r < i) 
(r > i) 


(19) 


When the system size N is odd, we diagonalize the Hamiltonian in the same way and 
obtain the energy gap as, 


A(iV,r,7) 


53 e ( fc s) - ^2 e ( k e) 

k$ ke 

2 ( r -!)+53 e ( k ^ - 53 e ^ 6 ) 

fc 5 k 6 


(r < i) 
(r > i) 


( 20 ) 
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where 


, „ , 2 . JV-3 , N-l 

h = 0, ±-tt, ■ • • , ±-^r, ±-^, 

, 1 3 N — 2 

k ‘ = ^ ^' " ’ ± m~ n ’ 


( 21 ) 

( 22 ) 


are wave numbers in the subspaces with odd and even parities, respectively. 

III. ENERGY GAP AS A FUNCTION OF THE SYSTEM SIZE 

In this section, we derive the asymptotic form of the energy gap and show that it behaves 
in an irregular way as a function of the system size. 


A. Isotropic case (7 = 0) 

We start with the isotropic case of 7 = 0, for which the Hamiltonian is already diagonal 
in Eq. (J5]), 


H — Y \^-2a\a k (cosk + T) + V 


(23) 


where the values of k depend on the boundary condition as in Eqs. ([6]) and ((7(). We find it 
convenient not to use the general expression derived in Eq. (fl9|i here, but we start with the 
above formula instead. For simplicity, we only consider even N here. 

We carefully identify the ground state in preparation to evaluate the energy gap. An 
apparent candidate for the ground state is the state where all modes with wave numbers 
satisfying 


cos k > —r, 


(24) 


are occupied. 

We check that those states have the number of fermions N c matching the respective parity 
sectors. Indeed, as should be clear from Fig. [2j the series of wave numbers k\ in Eq. ([6]) 
is consistent with odd N c and &2 is consistent with even N c . Notice that the parity of N c 
coincides with the parity of S^ ot up to a constant N/ 2, 

N N / . 1 \ N 

(25) 




3 = 1 


3 = 1 


7 






which is related to the statement that we consider the gap between states with even and 
odd values of S£ ot . 



FIG. 2: Left panel: Wave numbers for the presumed ground state in the case of k\ on the 
complex-fc plane. The black points mark the occupied modes, and the white ones are for 
the unoccupied ones. The number of occupied modes is odd. Right panel: Likewise, the 

number of occupied states for &2 is even. 


1. Energy gap 


The lowest energy in each parity sector is given by Eq. 
Eq. (j24li and a\ak = 0 otherwise. The energy gap reads 


, with a' k cik = 1 for k satisfying 


A(iv,r,o) = 


[—2 (cos k + T)] — [—2 (cos k + T)] 




(26) 


where ^7r and N - 7r are, respectively, the largest k\ and /c 2 satisfying Eq. (124)1 . Then, 


cos k = 1 + 2 cos -ji- 


|fcl|<^ 3=1 


N 


7T = 


(2 n + 1)tt 

Sin ^L. 

. 7T > 

Sill — 

N 


Sill 


cos k =- 


2mn 

~N~ 


\k2\<^ir 


Sill ■ 


7r • 

N 


(27) 


Next, we find n and m, as 


Nx 


Nx 1 


, m = 

- 1- 

[ 2tt J 


|_ 2tt 2 J 


n = 


x = arccos(—T). 


(28) 




















Dividing the quantities in the floor symbol above into integer and fractional parts, 


Nx 

77 

Nx 1 

77 + 2 


Nx 

27 


Nx 1 

27 + 2 


+ $i — n + hi (0 < hi < 1), 

+ S 2 — m + S 2 (0 < h 2 < 1). 


we find that for 0 < hi < |, 5 2 = hi + and the energy gap is 

r cos (^(1 - 4hi)) + 7l - r 2 sill (1 - 4hi)) 


A(iV,r,0) = 2 


cos ■ 


7r 

2N 


while for | < hi < 1, h 2 = hi — | and the gap reads 

r cos (^(3 - 4hi)) + 71-T 2 sin (3 - 4hj)) 


A(lV,r,0) = 2 


cos 


TT 

2N 


-T 


(29) 

(30) 


(31) 


(32) 


In both cases we can expand the above expressions in the powers of 1/A, obtaining the 
asymptotic formula 


A(A,r,0) 


7T 


vT^T 2 


N 


2 mod 


Ncj) 


IT 


,1 


+ 0(N~ 2 ), 


(33) 


where the angle 0 is defined as 


0 = arccos(r). (34) 

The standard, polynomial decay of the gap is modified here by (2 mod (^r, l) — l), which 
is a piecewise linear, continuos function of <f>N, oscillating between —1 and +1. 

In the Appendix A we present another derivation of the above formula using more general 
approach, suitable for anisotropic case as well. 


2. Anomaly in the energy gap 


We plot typical examples of the energy gap as a function of the system size in Fig. |31 
which have been calculated using the exact formulas in Eqs. (l3Tf and (l32j) . 

Most interesting, these figures look qualitatively very similar to those for the infinite-range 
quantum XY model [231] . 


N 


Aon = 


— Y( 

4A ^ v 


a ~ a j + 77/] 


i,j =1 


r N 

-Y 

2 ^ 

i =1 


<T, 


, N 

ll v "\ 
2 

i =1 


cq 


(35) 
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A A 



FIG. 3: System-size dependence of the energy gap for the isotropic XY model in Eq. (JT|) 
with 7 = 0 and h = 0. Left panel: Green, red and blue dots show values of the energy gap 
for T = 0.1, T = 0.3 and Y = 0.5, respectively. Right panel: The energy gap for Y = 0.999. 


For convenience, we reproduce some of them in Fig. [4j 

It is worth remaining that, in the infinite-range XY model, the gap can be tuned to close 
oolynomially, exponentially, or factorially by using an appropriate series of system sizes 
23 )]. Comparison of Figs. [3] and 0 strongly suggests that the same may be true for the 
one-dimensional model as well. We discuss it in some detail below. 


A 



FIG. 4: System-size dependence of the energy gap for the infinite-range XY model in Eq. 
(T35l) for the longitudinal field h =0 and various, fixed values of the transverse held T. Red, 
blue and green dots show the data for Y = 7r/10, Y = 1/3 and T = 1 — 7r/300, respectively. 
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3. Sensitivity of the energy gap to the change of magnetic field T 


There are several important points to notice in these results. Firstly, a slight change in 
T can lead to drastically different behavior. In order to illustrate this, it is useful to reduce 
the exact formulas for the gap in Eqs. (13TT) and (1321) for some specific values of T. 

For example, for T = x = arccos(—T) = 27 t/ 3, and hence Nx/2tt = N/ 3. Then, for 
N = 31 (l G N), hi = 0 according to Eq. (129]) and the gap takes a simple form 


A(IV = 3Z,T = 1/2,0) = y/3 tan (36) 

On the other hand, when IV = 3/ + 1, hi = 1/3 and for N = 31 + 2, hi = 2/3 and the gap 
becomes respectively 


A(7V = 3Z + l,r = 1/2,0) 


A (IV = 3/ + 2,r = 1/2,0) 


1 7 r \/3 7T 

- cos-sm- 

2 6IV 2 6 N 



1 

2 


1 

2 


(37) 


(38) 


For large N, Eqs. (I37H38P have the same system-size dependence since cos(7t/6IV) 
sin (7t/6IV). Thus, the gap follows essentially two separate curves as a function of the system 
size, given by Eqs. (135|) and (157|) (or (135]) ). as can indeed be seen in Fig. [3] for T = 0.5. 

A similar argument applies to x — arccos(—T) = nl/j with l,j G N satisfying 1/2 < 
1/j < 1, the latter condition coming from 0 < T < 1. In such a case 


Nx _ IN 

27“ 2j’ 


(39) 


and hi assumes a fixed value for selected series of system sizes. For instance, for N — ji (i e 
N), hi = 0 or |, and for N = ji + 1, hi = l/(2j) or 1/2 + l/(2j). Each of the series N = ji 
and N = ji + 1 gives a smooth curve describing A. Similar argument applies to other series 
of the form N = ji + const. Therefore, for some rational values of arccos(—r)/7T, the gap, 
as a function of the system size, follows a finite number of different curves. For irrational 
arccos(—r)/7T this argument does not apply and the gap behaves less regularly. 
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4- Envelope of the energy gap 


Secondly, the envelope of the gap as a function of the system size is proportional to its 
inverse. This is apparent when looking at the asymptotic expression in Eq. (]33lh and we 
show it as well in Fig. 0by plotting TV A as a function of TV for T = 0.1. 



FIG. 5: System-size dependence of the gap multiplied by the size for 7 = 0 and T = 0.1. 


This power-law dependence is indeed what is expected from the second order phase 
transition. One can easily construct the series of the system sizes which exhibits such a 
behavior by sticking to TV’s for which the oscillating term in Eq. (1331) is of the order of one. 
This would be a typical behavior as well. 


5. Possibility of rapid decay of the gap 


Note however, that we may be able to choose such a series of system sizes that the gap 
closes very rapidly. The reason is that the expressions in Eqs. (1311) and (j32j) oscillate between 
positive and negative values as illustrated in Figs. [3] and in the asymptotic formula in Eq. 
(l33]h This implies that the gap is crossing zero if TV is allowed to take continuous values. 
Then, it may be possible to choose appropriate values of TV close to those zeros which give 
very small, possibly exponentially small, values of the gap. If this is correct, which we failed 
to prove rigorously (unlike the case of the infinite-range quantum XT' model 23j), the rate 
of the gap closing could be tuned to behave exponentially by an appropriate choice of the 
series of system size, or be outright zero for some specific system size and magnetic field T. 
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6. Special case: independence of the system size 


Finally, Eq. (125]) simplifies considerably when cos(7r /N) < T. In that a case, the sum¬ 
mations over k\ and k -2 in Eq. (125]) run over all allowed values of wave numbers except for 
ki = tt and the gap reduces to a simple form l34j . 


A(iv,r,o) 


[-2 (cos k + r)] - 2(1 - r) - ^ [-2 ( cos k + r )] 

fcl k2 


2(1-r). (40) 


This last expression is independent of N as can be observed in the right panel of Fig. [3] in 
the range of cos(7t/A) < T. 


B. Anisotropic case (0 < 7 < 1) 


In the anisotropic case, we derive asymptotic formulas for the energy gap and show that 
its behavior differs remarkably in various regions in the ( 7 ,T) plane. Most intriguingly, we 
observe that the energy gap is oscillating in the incommensurate ferromagnetic phase, which 
is neighboring the isotropic critical line discussed earlier. We are able to link the period of 
those oscillations with the oscillations of the correlation functions in those cases. 


To that end, we start with Eqs. (1I9|) and (120]) and extend the procedure used in Ref. [261 
for the case of the Ising model, see Appendix A for details of the derivation. We summarize 
the results below. 

In the the ferromagnetic phase for T <1, which is the most interesting for us, 


A 


1/2 


-s 7 (r^ + 7 -1 - r 7 ypT7^I) f ' V* (1 + 0(Arl))> 


7r(l — y 2 ) 


0, 


4y/2 / 7 2 (l-r 2 )(l-r 2 - 7 2 ) 


7 r \ 


(I-7 2 ) 


1/4 j\r 

a 


Vn 


Vn 


cos(fjN + -0 o /2) (l + 0(N~ 1 )) 


(r 2 + 7 2 > 1 ), 

(F 2 + 7 2 = l), 

(r 2 + 7 2 < 1 ). 

(41) 


In the paramagnetic phase for T > 1, 

'87 (r 2 + 7 2 -1 - 7 r^r 2 + 7 2 -i) \ 1/2 A 


A w 2(r - 1) + 


7r(l — 7 2 ) 


Vn 


(1 + 0(N~ 1 )) 1 (42) 
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and, finally, on the critical line for Y = 1, 


A 


77T 

2 N 


2 7 - 


7 r 


+ 0(N~ 5 ). 


27 y 481V 3 

For convenience, a, A 2 , i/j and appearing above are defined as in Ref. 


(43) 


3fl] 


a = 


cos^ = 


1~7 

I + 7 ’ 

r 


(1 - 7 2) 2 

^0 = arg(y/r^ 


Ao — 


r 2 - 7 2 - 
r- [r 2 + 7 2 -i]^ 

1-7 




(44) 

(45) 

(46) 

(47) 


From these results, we notice that the energy gap oscillates as a function of N for T 2 + 7 2 < 1, 
which means that the ground state changes between subspaces with even and odd parity of 
N c as the system size N is increasing. On the other hand, the ground state for T 2 + 7 2 > 1 
always has even parity. (For Y > 0, 7 > 0. For negative T or 7 one can easily map the 
model onto the one with positive values of the parameters by suitable spin rotation. This 
leads to change of parity of the ground state for odd N and Y < 0). 

We plot the typical examples of the energy gap as a function of the system size in the 
ferromagnetic phase in Fig. [ 6 ] which have been obtained by numerically evaluating Eq. (1T9l) . 


A 



FIG. 6 : System-size dependence of the energy gap for an anisotropic XY model in Eq. (JT]) 
for h = 0 and in the ferromagnetic phase. Blue and red dots show the energy gap for 
(r, 7 ) = (0.85,0.95) and (r, 7 ) = (0.2,0.15), respectively. In the second case the gap falls 
toward zero very rapidly at periodically observed dips. 
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For (r, 7 ) = (0.85,0.95), i.e. r 2 + q 2 > 1, the gap decays in a standard, exponential way. 
By contrast, for (r,q) = (0.2,0.15), i.e. r 2 + q 2 < 1, there are additional dips visible, which 
mark the oscillations. 


IV. RELATION BETWEEN ENERGY GAP AND CORRELATION FUNCTIONS 


We notice that the connected correlation function in the infinite system, 


Of = (<K +R ) - (K)w 


i+H ), 


(48) 


has similar qualitative behavior as the gap. We focus on the ZZ correlation function here, but 


this observation is independent of the specific choice of local observab 


we quote the asymptotic of the correlation function, following Ref. 
In the ferromagnetic phase, T < 1, 

1 


30]. 


es. For convenience, 


2vr 


A - 2 *- 2 (i + OiR- 1 )) , 


0, 


(r 2 + 7 2 > i), 

(F 2 + 7 2 = l), 


Gf « < 


4 a 2R 

— -—SR 

71 R 2 


a iip(R+ 1 ) 


3? 


Jip(R- 1 ) 


1 - e 2 ^ \ 1/2 

1 — ct 2 e -2 *^ J 

1 - « 2 e- 2 ^\ 1/2 


_ g—2*7 


(l + ofi?' 1 )), (r 2 + 7 2 < i), 


where 3ft denotes the real part and a, A 2 and -0 are defined in Eqs. (fflll - fTTH) . 
In the paramagnetic phase, V > 1, 


G 


R 


if ( 1 + °(«- 1 ))' 


For the critical line, T = 1, 


G 


R 


7 r 2 R 2 

and finally, for the isotropic case of 7 = 0 and T < 1 

1 


(1 + 0 (R~ 2 )) , 


syR _ 

2 “ 7r 2 R 2 


sin(0R) 2 . 


(49) 


(50) 


(51) 


(52) 


Comparison of Eqs. (l4ll) - (l43D and (I49l) - (l52l) reveals close similarities. In the ferromagnetic 
phase, which is a first order transition from the viewpoint of our quantum annealing, the 
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gap disappears in an exponential way on the length scale given by the correlation length in 


the system, which, of course, agrees with the standard qualitative prediction 


25 1 . 


Note however, that the period of oscillations of the gap as a function of the system size 
in the incommensurate phase (y 2 + T 2 < 1 ) coincides precisely with period of oscillations of 
the correlation function. Remarkably, the same observation holds true for the second order 
transition in the isotropic case. Indeed, for 7 = 0, 0 in Eq. (133|) is exactly equal to 0 in Eq. 
(l52]h Trivially, this relation is also satisfied in other discussed cases, where neither the gap 
nor the correlation functions oscillate. 

For the XY model discussed in this article this relation is not accidental and can be 
understood since both the correlation function and the energy gap in the finite system are 
closely related to the Fourier coefficients of cos 9 k (sin 6 k ) in Eqs. (1 1 OH 11 [) and of e k in Eq. (fT2lh 


respectively (see 


30] and the Appendix A). Asymptotic behavior of the Fourier coefficients 


in both cases is determined by the poles of e{k) in the complex plane, namely by A 2 . When 
A 2 is real and positive there is no oscillatory behavior. On the other hand, for complex A 2 , 
its absolute value, a, is giving the rate of decay and its phase, 0, determines the period of 
oscillation. 

It is an interesting question if such a relation holds true also for other models, especially 
since it can severally spoil the performance of the quantum annealing protocol. While not 
directly relevant for us, we note, that for the systems in the thermodynamic limit there is 
a general and similarly looking relation linking the minima of the dispersion relation of the 
translationally invariant Hamiltonian (i.e. its spectrum) to the oscillations of the correlation 


function in the ground state 


39|. 


We discuss some further similarities between the gap and the correlation function in the 
finite system in the Appendix B. In particular, if T 2 + y 2 < 1 and the Fermion number is 
odd, a direct relation holds between the gap and the correlation function for finite size N 
as follows, 


of(Jv,r, 7 ) = (4K +R ) - <7)« +R ). (53) 

ap2N, r, 7 ) = (Wk. A (jv, r, 7)] . (54) 

This establishes a direct nontrivial relation between the correlation function and the energy 
gap. 
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V. CONCLUSION 


In this papers, we have analytically studied the energy gap of the one-dimensional quan¬ 
tum XY model with s = 1/2. 


We first analyzed the energy gap of the isotropic XY model to show that the energy 
gap at the second order transition point behaves quite anomalously as a function of the 
system size toward the thermodynamic limit. This resembles the property of the energy gap 
of the infinite-range quantum XY model, where the gap decreases in many different ways: 
polynomially, exponentially, or factorially, depending on the choice of the series of system 


sizes [23|. Also, our expressions of the gap, Eqs. 
conjectured from numerical calculations in Ref. 


(13TT) and ()32lh reproduce some of the results 
35] . In particular, their Eqs. (B3) and (B4) 


are corroborated by the asymptotic expansion in Eq. 


It is an interesting problem to study what aspects of these one-dimensional and infinite- 
range models are the key to the anomalous size dependence of the gap. The continuous 
symmetry is one of the common features of these models, but we should carefully investigate 
if this is the most essential property. From the viewpoint of implementation of quantum 
annealing, it is reassuring that the energy gap can be tuned to decrease polynomially as a 
function of the system size even when the lower bound closes very quickly, because the right 
choice of the series of system size allows us to avoid the problematic rapid closing of the 
gap. This also means that that one should be careful to interpret the results obtained from 
numerical simulations or other methods for a limited series of system sizes because different 
properties may emerge for different series of system sizes. 


Finally, we have analyzed the energy gap of the anisotropic XY model at the first order 
transition point. We obtained closed asymptotic expressions for the gap and notice sim¬ 
ilarities between the energy gap and the correlation function. For instance, both for the 
anisotropic and isotropic case the oscillations of the energy gap in the finite system precisely 
coincides with the oscillations of the correlation function in the thermodynamic limit. It is 
an interesting future problem to investigate whether or not similar relations hold in other 
models. 
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Appendix A: Evaluation of the energy gap 


For the Ising model (7 = 1), a closed analytical form of the energy gap is obtained in 


Ref. 


26]. We generalize the method to the XY model. 


Expanding e(k) appearing in Eqs. (f2Ul) and (ITTTj) as a Fourier series 0,1 cos(Zfc), we 
can write the energy gap as 26]: 


A = 


-2N (a/v + a 3 N + ®5N + •••)! 


(r < 1 ) 


2 (r — 1 ) — 2 iv(ajv + CL3N + « 5 n + • • •) (r > 1 ) 


(Al) 


where 


az = - 
7 r 


dk cos (hi) yJ{T- cos k) 2 + y 2 sin 2 k. 


(A2) 


1. Ferromagnetic phase, T < 1, with r 2 + y 2 > 1 


If we introduce z = e tk , Eq. (IA2|) is equivalent to 


ai = l i I ( Z ~ Xl ^ Z ~ X ^ Z ~ Al ~ Aa 

^ J\z\=i i 

L-l \-l 


Az 2 


(A3) 


with branch cuts along (A x , A 2 ) U (A 2 , Ai), where Ai and A 2 are defined as follows, 


Ai = 
A 2 = 


r + [r 2 + 7 2 -i]^ 

1-7 

r- [r 2 + 7 2 -i]^ 

1-7 


(A4) 

(AS) 


By deforming the integration contour around the branch cuts, we obtain 


ai — 


y/T" 


7/■' 


7 r 


A 2 l I t l 2 (1 — t) 1//2 (t — A 2 /Ai) 1//2 (A 1 A 2 — t) 1 A(i — t/\iy/ 2 dt. (A6) 

J A 2 /A 1 
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For large Z, the integral is dominated by the contributions around t = 1. We further assume 
that we are away from the line T 2 + 7 2 = 1, that is A 2 /A x < 1, and expand the integrand 
around t = 1 taking the leading order in t — 1. We arriving at 


ai 


y/a 


T) ,-i, 


-A/(l - A 2 /Ai) 1/2 (AiA 2 - 1) 1/2 (1 - 1/A 2 ) 1/2 / t*" 2 (l - t) 1/2 dt 

71 J a 2 /Ai 

\-i ( 87 fr 2 + 7 2 -1 - 7 ryr 2 + 7 2 -1) \ 1/2 r i 
v --il+- ! ™ 


7 r 


We can evaluate the above integral as follows, 

[ t l ~ 2 (l-t) 1/2 dt& f t l ~ 2 (l-t) 1/2 dttt 

J \\/\2 j o 2 ly/l 

Therefore, 


dvr 


a -jv / 8 7 (r 2 + 7 2 -1 - 7 iyr 2 + 7 2 - 1 ) 


1/2 


7T 


and for N large enough, 


(I-7 2 ) 


(in ^ <25jv • • • 


2NVN 


Using the above results, we obtain the asymptotic behavior of the energy gap as 

87 fr 2 + 7 2 - 1 - 7 r^r 2 + 7 2 -i) \ 1/2 x -n 
A « —2 Ncln ~ 1 ^ / 1 2 


7r(l — 7 2 ) 


Vtv' 


(AS) 


(A9) 


(A10) 


(AH) 


2. Critical line, F = 1 


Taking the limit of T —> 1 in Eq. 

\/i - 7 2 r 1 


gives 


a* = 


(1 -‘> h i + 7 

1+7 


1- 7 A 1/2 / I + 7 


— t 


1/2 


dt. 


(A12) 


7T / 1-7 ' \ i + 7/ \ 1 — 7 

For large Z, the integral is dominated by the contributions around t — 1. We expand the 
integrand around t — 1 and take the leading order in f — 1, arriving at 


ai — 




a r 


7T 


'0 


27(1-*) 7(1 ~ ZQ 2 (1-t) 3 

\/l - 7 2 \/l - 7 2 4 7 a/1 - 7 2 


+C?((1 - Z) 5 )] dt + 
27 1 


7 2 V 1 - 7 2 

1_7 ^ ^(z- 1 ) 


87 


7 


1 + 7 
2 7 3 


7 r Z 2 — Z 7r(Z 3 — Z) 27 t 7 Z 4 + 2Z 3 — Z 2 — 21 717 Z 5 + 5Z 4 + 5Z 3 — 5Z 2 — 6Z 

7T Z 2 \27T7 7T y Z 4 1 


C>(Z- 6 ) 

(A13) 
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Therefore, we obtain 


A 


-2N J2 

s=1,3,5,... L 


2 7 1 
7 r s 2 N 2 


+ 


27T7 


27 
7 r 


s 4 iV 4 


C>(A 


- 6 \ 


77T / 3 

= 2iV + ( 7 ~ ^7 


7 r 


48A 3 


+ e>(A 


-5\ 


(A14) 


3. Paramagnetic phase, T > 1 


If we introduce z = e lk , Eq. (1A2|) is equivalent to 


a , = - r (z ~ Al)(z ~ A2)(z 2 ~ Ar ‘ )(g ~ . 
ft J\z\=l * ' 

with branch cuts along (A^ 1 , A 2 ) U (A^" 1 , Ai). 

By deforming the integration contour, we are lead to 


(A15) 


a, = -AZZv r' 


2 1 2 (1 — t) 1 ^ 2 (t — 1 /AiA 2 ) 1 ^ 2 (Ai/A 2 — t) 1 ^ 2 (l — A2t) 1 ^ 2 dt(Al6) 

77 J 1 /A 1 A 2 


For large Z, the integral is dominated by the contributions around t — 1. We expand the 
integrand around t = 1 and take the leading order in t — 1, arriving at 


a, « -A- 2 A i (1 _ l/X^y/^XJX, - 1 ) 1 / 2 (1 _ a=) 1/2 [' ( 1 - 2(1 _ 7 / 2 *; 

^ A1/A1A2 


A* /s 7 (r 2 + 7 2 - 1 - qTv/r 2 + 7 2 - 1) 


1/2 


7T 


1 — 7 2 

A ; /8 7 (r 2 + 7 2 -1 - 7 r ^r 2 + 7 2 - 1 )' 


'l/AiA2 


t /_2 (l - t) 1/2 dt 


1/2 


7T 


(I-7 2 ) 


2Za/T 


(Air) 


Similarly as for the ferromagnetic case discussed above, for N large enough the asymptotic 
of the gap is determined by a at. From the above results, we obtain the asymptotic behavior 
of the energy gap as 


A « 2(r- 1) + 



-1 - 7 ryr 2 + 7 2 -i) 

7r(l — q 2 ) 


1/2 


A? 

Vn' 


(Ais) 
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4. Incommensurate ferromagnetic phase, r 2 + 7 2 < i 


Using z = e zk , Eq. (IA2j) i 


is 


ai = 


i^/l 


7 


2 vr 


dz t (z — a;e ^) 1//2 [z — a. 1 e*^’) 1 ^ 2 (^ —ae U ’) 1 / 2 (z — a 1 e 


z 1 / 2 


;l /2 


r UI=i z 

with the branch cuts along ( 0 , ae l ^) U (o _ 1 e*0 oo) and ae 1 ^ — (T + i\J 1 — T 2 — 7 2 )/(l + 7 ). 
That is, the phases are dehned in the region (0, 27T + 0) for the branch cuts in the upper 
half-plane, and (— 0 , 2tt — - 0 ) for the cuts in the lower half-plane. 

By deforming the integration contour, we find 


ai = 3ft 


2 iy/l — 7 2 f 1 dt 
f 2 


( ae l ^t ) (1 — t) 


- /W 2 


7T 7 0 

(a" 1 - ot) 1/2 (te^ - e -^) 1/2 (toe^ - cTV^) 172 ' 


(A 20 ) 


To estimate the leading behavior in N, we expand the above integral around t — 1 to get 


a N pz 


2i ^ 7 V e ^ f 1 P N (1 - t ) 1 / 2 («' 1 - a ) 1 / 2 ( e * - e ^) 1/2 1/2 


ae - a e 


*)’ 


= ^ a w cos(ipN + */ 2 ) [ *(1 - t )VV 

4V2 /V(l - r 2 )(l - r 2 - 7 2 )\ 1/4 a" ,. , 7 , ,. 

Therefore the asymptotic of the gap for large enough N, 


N -2 


(A 21 ) 


A 


4-02 07 2 (1 — T 2 )(l — T 2 — 7 2 )\ 1/4 a N . , Ar , 

— (--W--) + W 2) 


(A22) 


5. Isotropic critical line 7 = 0 


To that end we can take the limit 7 —> 0 in Eq. (1A20I) . which is equivalent to setting 
a = l and 0 = 0 = arccos(T). 




a/ = 5R 

This integral can be calculated and subsequently expanded for large l as, 
ai = 3ft 


(A23) 


4i(cos0-flsin0) l(j)l 
7 r(Z 2 - 1)/ 


/1 r 2 r 

~ -4--— cos(/0) + 4— sin (70) + C>(Z -4 ) 

7U Z 7U' 3 


(A24) 
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Now, in the leading order in N the gap is 


A ~ 


yi - r 2 ^ cos(siV0) 

2-^ v 


71N 


(A25) 


s=l,3,5,... 

The sum in the above expression gives the oscillating term periodic in Ncj) with the period 
27r. Calculating the sum, we obtain 

nVl - r 2 ( 2 


A — 


N 


IT 


arcsin(cos(Ad)) 


(A26) 


which can be also expressed using modulus as in Eq. 

To complete the analysis we consider the next term in the expansion in powers of 1/A. 
Extending Eq. (1A25|) we get 


A ~ 


, Vl - r 2 cos(sA0) 

2 A v ^ 


7 TN 


7 tN 2 


E 


sin (sN(f>) 


(A27) 


s=l, 3,5,... s=l,3,5,... 

The second term is also periodic in Ncj) with the period 27 t, but the zeros are shifted by 
7t/2 . It is piecewise quadratic in <f>N. The second term may also become more important 
when approaching the critical point at T = 1, as suggested by the prefactors. Performing 
the sums, and writing the result in a compact form we obtain 

tiVI - r 2 ( 2 


N 


IT 


arcsin(cos (Ncj))) 


tt'V f 2 / 2 

H—— ( — arccos(cos(A</>)) ( 2-arccos(cos(A</>)) Sign(sin(A0)) 

AN 2, V 7T V 7T 


For instance, this shows that the expansion in Eq. 
of zeros of the energy gap. 


(A28) 

does not capture the exact position 


Appendix B: Relation between the gap and correlation function in the finite system 


We consider the ZZ correlation function of the one-dimensional quantum XY model in 
a finite system in the ground state. It has been calculated in Ref. |30|, 


G?{N,T,i) = 

= [— sin 9k sin (Rk) + cos 9 k cos (Rk)] 

k 

x [sin 9k sin (Rk) + cos 9k cos (Rk )], 


(Bl) 
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where 


and 



2 N - 2 

0 ' ± N 7! ’'"’ ± ^r 7! ’ T 

1 3 TV — 1 

± N 7r ' ± N ,T -'-- ± — 7 ' 


(N c odd) 
(N c even) 


(B2) 


cos 6 k 
sin 0 k 


cos k + r 

-^/(cos k + T ) 2 + (ysinfc ) 2 
7 sin k 

■\J (cos k T ) 2 + (ysinfc ) 2 


(B3) 

(B4) 


Let us first consider the region r 2 + y 2 < 1. We now show that the correlation function 


G ??(2 N, r, 7) is written explicitly in terms of the energy gap if T 2 + y 2 < 1 and the Fermion 
number is odd. In this case, because shi^Nk^) = 0 in Eq. (IBID , the correlation function 
simplifies to 


Gf(2JV.r, 7 ) 



(B 5 ) 


where 


k 7 = 
h = 


1 2 


N - 1 


N 77:77 


2 N — 2 

0, ±-7T, • • • , ±-7T, 7T 

N ’ N 


k2 


1 3 

± iv 7r ’ ± ]v 7r ’ 


N -l 


N 


-IT. 


(B6) 

(B 7 ) 

(B8) 


Furthermore, according to Eqs. (fT2l) and (1B3jl . we find cos 9 k = de{k)/dT, and the correlation 
function is seen to have the following expression, 

G72/v.r, 7 ) = J i 5 ^^ £ (M-E f (fe))) • ( B9 > 

In the region 0 < Y < 1, from Eq. (|T9lh the gap is 


We therefore have 


A(iV,r,7) 


fcl fc 2 


Gf(2fV,r,7) 


/ 1 d 

V2iV9f 


A(iv, r, 7 ) 


2 


(BIO) 


(Bll) 
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It is impossible to find a similar relation for an even number of Fermions because sin (Rk) 
in Eq. (IBID does not vanish. Likewise, in the region T 2 + 7 2 > 1, because an even number 
of Fermions always gives the ground state, we are not able to find a similar relation. 
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